Steep sharp-crested gravity waves on deep water 
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A new type of steady steep two-dimensional irrotational symmetric periodic gravity waves on 
inviscid incompressible fluid of infinite depth is revealed. We demonstrate that these waves have 
sharper crests in comparison with the Stokes waves of the same wavelength and steepness. The 
speed of a fluid particle at the crest of new waves is greater than their phase speed. 

PACS numbers: 47.35.+i 



A proper understanding of various wave phenomena on 
the ocean surface, such as modulation effects and insta- 
bilities of large amplitude wave trains [Q , , formation 
of solitary freak and breaking waves |^, re- 
quires knowledge of a form and dynamics of steep water 
waves. For the first time surface waves of finite ampli- 
tude were considered by Stokes 0. Stokes conjectured 
that such waves must have a maximal amplitude (the 
limiting wave) and suggested that a free surface of the 
limiting wave near the crest forms a sharp corner with 
the 120° internal angle (the Stokes corner flow). A strict 
mathematical proof of the existence of small amplitude 
Stokes waves was given by Nekrasov. Toland Q proved 
that Nekrasov's equation has a limiting solution describ- 
ing a progressive periodic wave train which is such that 
the flow speed at the crest equals to the train phase speed, 
in a reference frame where fluid is motionless at infinite 
depth. Longuet-Higgins and Fox constructed asymp- 
totic expansions for waves close to the 120°-cusped wave 
(almost highest waves) and showed that the wave profile 
oscillates infinitely as the limiting wave is approached. 
Later, in the crest of a steep, irrotational gravity 
wave was theoretically shown to be unstable. 

The purpose of the present work is to give evidence 
that a second branch of two-dimensional irrotational 
symmetric periodic gravity waves of permanent form ex- 
ists besides the Stokes waves of the same wavelength. 
The original motivation is as follows: the Bernoulli equa- 
tion is quadratic in velocity and admits two values of the 
particle speed at the crest. The first one corresponds 
to the Stokes branch of symmetric waves for which the 
particle speed at the crest is smaller than the wave phase 
speed. The opposite inequality takes place for the second 
branch which might correspond to a new type of waves. 
In the second part of the Letter, we prove this numeri- 
cally by using two different methods. 

Consider a symmetric two-dimensional periodic train 
of waves which propagates without changing a form from 
left to right along the x-axis with the constant speed c 
relative to the motionless fluid at infinite depth. The set 
of equations governing steady potential gravity waves on 



a surface of irrotational, inviscid, incompressible fiuid is 
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-oo < y < 77(6*); (1) 
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9 = x~ct. (4) 



where $(6', y) is the velocity potential, r]{9) is the eleva- 
tion of a free surface, and y is the upward vertical axis 
such that y = is the still water level. We have chosen 
the units of time and length such that the acceleration 
due to gravity and wavenumber are equal to unity. 

As it follows from the Bernoulli equation (||) , a solution 
may be not single-valued in the vicinity of the limiting 
point. Indeed, the particle speed at the crest q{0) is 
horizontal and is defined as follows: 



(0, 7j{Q)) = q{0) = c ± - 277(0), 



(5) 



77(0) being the height of the crest above the still water 
level. The "— " sign corresponds to the classical Stokes 
branch. The value 77inax(0) = c^/2 corresponds to the 
Stokes wave of limiting amplitude. In this case, the par- 
ticle speed at the crest is exactly equal to the wave phase 
speed: (7max(0) — c. Taking into account both signs 
in expression (||), we assume that a second branch of 
solutions should exist apart from the Stokes waves, at 
77(0) < 7;inax(0). The particle speed at the crest of a new 
gravity wave must be greater than c and has to increase 
from c to 2c while the wave height decreases from 77inax(0) 
to 0. Moreover, the mean levels of these two fiows relative 
to the level y = of still water must also be different: 
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Thus, the existence of a second branch of solutions of the 
set of equations (|l|)-(§) does not contradict Garabedian's 
theorem that gravity waves are unique if all crests 
and all troughs are of the same height because the latter 
was proved for a flow with the same mean level. 

To construct a numerical algorithm we use the method 
of the truncated Fourier series and the collocation 
method, in a plane of independent spatial variables. 

The method of the Fourier approximations. Let us in- 
troduce the complex function R{9, y) such that 
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-ic{R-R*), '^ = c{R + R*) 



(7) 
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where ^' is the stream function, * is the complex con- 
jugate. Using the relations ~ = —'^x, the 
kinematic boundary condition can be presented as 
follows: 



dx 



{R (9, T]) + R* (61, r/) - rj{0)) = 0. 



(8) 



Approximate symmetric stationary solutions of Eq. (0, 
H, H, ^ are looked for in the form of the truncated Fourier 
series with real coefficients 



N 



R{e, y) = ^enexp(n(2/ + z0)); 



71 = 1 

M 



X! '7nexp( 



mt)), r]-n = rin, 



(9) 



(10) 



71 = -M 



where the Fourier harmonics ^„ , , and the wave speed 
c are functions of the wave steepness A determined by 
the peak-to-trough height: 



: = - > ?72n+l, 11 

n=0 



2tt 



square brackets designate the integer part. Substitution 
of expansions (^ and ( px| ) into the dynamical and kine- 
matic boundary conditions (^, (0) (the Laplace equa- 
tion (^ and boundary condition (^ are satisfied exactly) 
yields the set of -I- Af + 1 non-linear algebraic equations 
for the harmonics 77„, and the wave speed c 
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n2 = l 



where f^^ are the Fourier harmonics of the exponential 
functions exp{nir](9)): 
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exp(ni7?(0)-*n0)d0, ^ f^K (14) 



They were being calculated using the fast Fourier trans- 
form (FFT). In addition to these equations, the connec- 
tion (|ll|) between the harmonics i]n and the wave steep- 
ness A should be taken into account. 

The set of equations (|l|), (|l|) was being solved 
by Newton's iterations in arbitrary precision computer 
arithmetic. Since the non-linearity over ^„ and rjn is of 
a different character (polynomial and exponential), the 
value of M should be chosen greater than N to achieve 
good convergence. A different number of modes for the 
truncation of the Fourier series (||) , ( p^ was also used by 
Zufiria Ifl2| in the framework of Hamiltonian formalism. 



The method of collocations. The harmonics ^„ of ex- 
pansion (^ can also be found in another way without 
expanding elevation into the Fourier series. In this ap- 
proach, Eq. (H) and explicitly integrated Eq. (H) are 
to be satisfied in a number of collocation points 6j — 
jn/N, j — 0, N, equally spaced over the half of one 
wavelength from the wave crest to the trough, similar to 
Rienecker and Fenton |l|]. This leads to 2 A^ -I- 2 alge- 
braic equations for the harmonics the values of the 
elevation 77 at the collocation points, and the wave speed 
c. To make the numerical scheme better convergent, the 
greater number of collocation points may be used in the 
dynamical boundary condition (^: M = PN, P is an 
integer. 

The results of calculations and discussion. The de- 
pendence c{A) of the speed of steep gravity waves on 
their steepness is shown in Fig. Along with the curves 
obtained by the Fourier and collocation methods, we in- 
cluded high accuracy calculations of the Stokes branch 
by the method of an inverse plane according to the equa- 
tions presented in Tanaka's paper |lj]. In the plot, 
point I (A — 0.13875) is the maximum of wave speed, 
point 2 {A = 0.14092) is the relative minimum, point 
3 {A — Ayna,x = 0.141074) corresponds to the limiting 
steepness at A^ and M given. For greater values of A^ 
and M, Amax > 0.14108 is obtained. Note, that less 
accurate calculations by the collocation method give a 
greater value of the limiting steepness which is close to 
that reported by Schwartz [|l5j . 

The key result of our numerical investigation is that we 
have revealed a new branch which arises from the point 
of the limiting steepness in the direction of its decreas- 
ing so that the loop 2-3-4 is formed. Thus, the point 
of the limiting steepness seems to be the point of maxi- 
mum of A, not the breaking point. It should be noted, 
that firstly we obtained the new branch by the Fourier 
method, and only after that we could track it by the col- 
location method using the starting points generated by 
the first method. As it is shown in Fig. |[ the profile of 
the new solution near the crest is sharper than the profile 
of the Stokes wave of the same steepness. Because of this 
we named it "the spike wave". The difference between 
the crests of the Stokes and spike waves becomes stronger 
as wave steepness drops relative to the limiting value. In 
Fig. ^ the dashed lines designate the exact local Stokes 
solution (the Stokes corner flow) which corresponds to 
the limiting wave with a maximal value of A (point 3 in 
Fig. 1^). In the immediate vicinity of the crest, the profiles 
of the almost highest Stokes and spike waves asymptoti- 
cally tend to the dashed lines. This tendency is seen to 
have the oscillatory character for both waves. For the 
Stokes waves such oscillations were analytically obtained 
earlier in 

The values of the wave speed c for the Stokes and spike 
waves calculated by the Fourier method at different val- 
ues of wave steepness are presented in Table ||. The 
high accuracy values of the wave speed for the Stokes 
branch obtained by using Tanaka's procedure are also in- 
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FIG. 1: The wave phase speed c of a surface wave versus its steepness A. 



TABLE L The values of the wave speed c and the mean water level r/o for the steep Stokes and spike waves with the wave 
steepness A calculated by the Fourier method. The values without footnotes were calculated at A*' = 200, M = 3N. The wave 
speed for the Stokes waves obtained by Tanaka's procedure are presented to estimate the accuracy of our calculations. 







Stokes wave 




spike wave 




A 


c 


C Tanaka 


r,(^' X 10^ 


c 


77^^) X 10^ 


0.14 


1.0926149034 


1.0926149034 


-2.46-10"^'' 


1.09246 


-130 


0.1406 


1.09233763 


1.0923377499 


-1.3M0"* 


1.09249 


-60 


0.14092 


1.0922742 


1.0922768392 


-0.0806 


1.092422 


-22.3 




1.0922761" 




-0.0221 


1.092427" 


-17.0 


0.141 


1.0922796 


1.0922808596 


-0.385 


1.092389 


-15.3 


0.14106 


1.0922949" 


1.0922851047 


-1.20 


1.0923548" 


-8.28 




1.0922962' 




-1.07 


1.0923550' 


-7.40 


0.14107 


1.0923008" 




-1.64 


1.0923458" 


-6.82 


0.14108 


1.0923114" 




-2.55 


1.0923321" 


-4.89 




1.0923145' 




-2.45 


1.0923303' 


-4.08 



"AT = 200, M = 47V. 
I'M = 210, M = 4N. 
"N = 250, M = 3N. 



eluded for comparison. One can see, that for the Stokes 
branch the accuracy of the Fourier method gradually de- 
creases as wave steepness increases up to the almost high- 
est steepness A = 0.14108. The correspondent value of 
the wave speed has only 5 digits stabilized. Note, that 



Tanaka's procedure diverges a,t A > 0.141062. While 
moving along the new branch the accuracy becomes still 
less, and much greater N are needed to stabilize a greater 
number of digits. As a result, the form of the loop in 
Fig. has not yet stabilized at = 200 and will en- 
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FIG. 2: The profiles of a free surface for the Stokes and spike 
waves of the same steepness and wavelength near the crest. 



TABLE II: The values of the particle speed at the crest of 
the Stokes and spike waves of the same wave steepness A, in 
a reference frame moving with wave speed. All values were 
calculated by the Fourier method at = 200, M = 4A. 



A 


9(0) - c 

Stokes wave 


spike wave 


0.14092 


-0.0370 


0.0540 


0.14106 


-0.0121 


0.0240 


0.14107 


-0.0076 


0.0194 


0.14108 


-0.0003 


0.0121 



large with increasing A^, the cross-section point with the 
Stokes branch being moved to the left. 

Table || also demonstrates that besides the form near 
the crest, the Stokes and spike waves of the same steep- 
ness have different mean water levels 770 relative to the 
still water level [see Eq. (||)]. One can see, that at 



the Stokes branch ?7q ' rapidly descends as A decreases, 
whereas rj^^^ increases for spike waves. Analysis of de- 
pendences of ryg^"* , TyQ^' on N and M indicates that they 
tend to different values at A^, M -^00. 

At the beginning of the paper we assumed the existence 
of a new type of gravity waves for which the speed of a 
particle at the crest is greater than wave speed. This 
property is confirmed by the calculations presented in 
Table |l| 

Thus, the spike waves, which we found numerically 
using two independent methods, present a new type of 
gravity waves we looked for. In the present work, we 
interested only in the existence of new stationary solu- 
tions and did not investigate their stability. Profiles of 
the almost highest Stokes and spike waves differ only in 
the vicinity of the crest. This leads us to an assumption 
that excitation of spike waves may possibly be connected 
with the crest instabilities |0 of the Stokes almost high- 
est waves. From the other side, sharpening of the crest of 
a spike wave, when wave steepness decreases (see Fig. H), 
makes us look for a relation to a problem of existence 
of solitary waves on deep water. At present, all existent 
experimental observations of surface solitary waves on 
deep water are usually interpreted by excitation of inter- 
nal waves in stratified ocean |p^ . However, verification 
of our assumption demands another numerical algorithm 
since the ones presented above become ineffective. Fi- 
nally, two- valued character of a solution of Eq. (|l|)-(0) in 
the vicinity of the limiting steepness does not depend on 
depth, as follows from Eq. (^) . We have recently revealed 
a second branch for a layer of finite depth. 
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